Tunneling of the Closed Friedmann Universe with 
Generation of Scalar Waves 



> 



X 



V. Ts. Gurovicht H.-J. SchmidtWd I. V. Tokareva * 



O ■ February 7, 2008 

O 
O 



Abstract 



The evolution of the closed Friedmann Universe with a packet of short scalar 
waves is considered with the help of the Wheeler-DeWitt equation. The packet 
ensures conservation of homogeneity and isotropy of the metric on average. It is 
shown that during tunneling the amplitudes of short waves of a scalar field can 
increase catastrophically promptly if their influence to the metric do not take into 
account. This effect is similar to the Rubakov-effect of catastrophic particle creation 
calculated already in 1984. 
| In our approach to the problem it is possible to consider self-consistent dynamics 

of the expansion of the Universe and amplification of short waves. It results in a 
decrease of the barrier and interruption of amplification of waves, and we get an 
exit of the wave function from the quantum to the classically available region. 
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1 Statement of the problem 

Here we present the basic equations of the problem. 

In the space-time of the closed cosmological Friedmann models with space-time metric 

ds 2 = dt 2 - a 2 (t)[d X 2 + sm 2 X (d6 2 + sin 2 6 dtp 2 )}, (1) 
where < Xi $ < 71 1 < <p < 2n , the action is assumed to be 

R - 1 , 1 



S 



A +-<p ,l <p,i - -m 2 (p 2 



16nG 2 r r ' 2 



-gdn, (2) 



Here R is a scalar curvature, ip is a potential of the scalar quantum field with a mass m, 

(H — c — 1), A is a cosmological A - term. 

The equation for ip has the form Dip + m 2 ip = 0. Below let us set the scalar field as 
radial standing waves with time-dependent amplitudes 

<Pk = Mt)f k (x); A(x) = — (3) 
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Note that the functions fk{x) do not have singularities on the poles (x = 0, ir) for integer 
numbers k, and that they are orthonormal on the coordinate x m the measure defined by 
metric (D), i.e., 

J sin 2 xfkfn dx = hn ■ (4) 

Now the action (Q) can be integrated over the angular variables, and we rewrite it 
after replacement of R by its value through a as 
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dt, (5) 



where 



A fc = ^8G/3A k (t); A = 8ttGA/3. 

These redefinitions are made to get G and tt only in front of the integral, and the 
integrand to be free of both constants. 

In this formula, the point denotes the derivative with respect to t, and one supposes 
that all k ^> 1, i.e., we assume that waves with wavelength comparable to the diameter 
of the Universe are negligible for the dynamics. 

The waves packet (J2k) ensures the homogeneity of the model during expansion with a 
high accuracy. After introduction of generalized momenta the Hamiltonian function can 
be presented as 

H = -^-a + Aa 3 + y(^ + k 2 A\a + m 2 a 3 Al \ = 0. (6) 
4a k\ Aa J 

Here P a , a are a generalized momentum and a coordinate for the Universe as entire; 
and Afc are the same for each radial mode of the field. Corresponding to eq. (^) the 
Wheeler-DeWitt equation for the wave function of the Universe (WF) has the form 

where the superpotential V is given as 
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Here we are neglecting J2k ^ 2 Af,a 4 in comparison with the sum residual in eq. (||) since 
k ^> 1. This is justified because in the short-wavelength approximation we are applying 
here, the mass of the particles does not have much influence. 

The factor ordering parameter p Q is also introduced in eq. (H). 

2 Semiclassical approximation 

We are going to search a solution of equation (|7]) in semiclassical approximation 

^ c = exp(*S c ), (9) 
where 5* c is a classical action. In this approximation we have 
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This nonlinear equation of the first order is similar to the Hamilton- Jacobi equation of 
analytical mechanics. For its solution one can take advantage of the complete integral 0] . 
It can be found from the system of characteristic equations [[| . As the eq. ( fLOf) supposes 
the solution for {OS/ da) , the scale factor becomes an argument along the characteristics. 
From here we have 

dS c _ _V_ dkk _ _ JZfc_. <kk _ 4fc 2 A fc a 2 
da F' da a 2 F' da F 

Here V is indicated by (|8|), and 

F = -^(£ql)l^-V. 

Now the evolution of the WF \1/ goes along a characteristic. Turning points are the 
roots of V(a) = 0. They separate classically available and forbidden regions along the 
characteristic eq. (|TTJ) Q . A similar approach for the strongly anisotropic WF in a Bianchi 
type I model is presented in ref. ||, and for a closed Friedmann model in ref. |J. 

In the region under the barrier the semiclassical solution can be found by an ansatz 
of the form 

* = exp(-|5 , e |). (12) 

The equation for S e can be obtained from eq. ( |T0D by changing V —V. Thus, as 
distinct from eq. (|TT|) the system of the characteristics has the form 

dS e _ V 
da F P ' 
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3 The evolution of the scalar field modes in the clas- 
sically allowed region 

Let us investigate how amplitudes of short wavelength change at the above mentioned 
evolution of the WF. Equations for q k and A k follow from ( [LTD 



da a 2 F' da 
where 



dA k q k dq k Ak 2 A k a 2 



F = -^(T / q 2 k )/a 2 -V. 

This system is an essentially non-linear one since co-factors of q k and A& in the right 
hand sides of the equations depend on a, too. However, the presence of the large parameter 
k 2 allows finding a solution. 

Following from the system (|T4|), the equation for the generalized momenta q k is 



d 2 q k Ak 2 q k (2 F'\ dq k 

da 2 F 2 [a F I da' 1 ' 



where ' = d/da. Here F(a) is an unknown function according to eq. (O). 



It is easy to show that the following solution 
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satisfies the last equation with the high accuracy of 0(l/k 2 ). And correspondingly, 
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sin $(a), 



(17) 



where Cko is an arbitrary constant, and the phase of the solution is chosen from the re- 
quirement of amplitude finiteness at a — > 0. The mentioned accuracy of solutions requires 
the assumption that the logarithmic derivative F'/F does not give a large parameter k. 



The last statement follows from the fact that each mode of wave enters into eq. (14) for 
t 1 as a term 
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It is clear from substitution eqs. fllq 
rewriting of the amplitude Cko/V^k 
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(0) in flU). With taking into account (|T 
Cfc, we have for F the following expression 
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Further, we suppose that A <C 1 and we have in the superpotential 
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Values of V and F define an action of WF \I/(a) according to the system of charac- 
teristics (|TTD. Because of the large number of short waves, the sum in eq. (|20| ) can be 
taken as fast oscillations average with an adequate accuracy. Then the expression for the 
superpotential gets the form 



V = 4 
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The existence of an inner classically available region requires a realization of the con- 
dition V(a) < in the mentioned interval < a < a . (Let us remember — see eq. (|6]) 
- that the total energy E vanishes for the closed universe.) The potential barrier should 
begin at a = a , i. e. V(a) > at a > a . Let us notice that the vanishing of the factor 
ordering parameter (p = 0) allows us to get an inner classically available region only at 
the expense of the sum J2 Cf • But in that case, a potential barrier necessary for a further 
development is absent. So, we choose the factor ordering p = 1 in the following. 

Owing to the amplitude of short waves begin to increase quickly under the barrier, 
we assume that their total energy is negligible in the considered region V(a) < 0, and 
the evolution of the Universe is determined only by the first two terms in the r.h.s. of 
eq. (pl|). In this case the dynamics of the amplitudes of the short waves can be given 
explicitly. The phase in eq. fllSD is 
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Here, ao can be obtained from the condition 



V{a ) = 4 
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A case when the total energy of the modes is not negligible in the expressions for 
V and F can be similarly considered. But then the final expression becomes essentially 
complicated. 



4 Modes of the scalar field under the barrier 

The superpotential V is positive at a > ao- This means that the evolution of the WF 
occurs under the barrier. At this range, the problem is defined by the system of charac- 



teristics (J13Q , and we have for q k the following equation 

da 2 F 2 [a FJ da' 



(24) 



This equation has an inverse sign at first term on the right in comparison with eq. 
(f5|) for the classically allowed region. It changes the solution in principle, but for a large 
parameter k 2 we can find a solution by a method similar to the one used for eq. ([15]) 

Qk = C k a cosh $ e (a); 

A fc = ^|sinh$ e (a); (25) 
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Here the conditions of continuity of q k and A^ at the point a = a were used for the 
selection of the constants of integration. Both in section 3 and here, the correctness of 
the solution ( ^5] ) is determined by the requirement of F' e /F e <C k. This requirement is 
checked by substituting of the obtained solution in F e of eq. (|13|). There 



Fe = -p2C 2 + Aa 2 . 

In the last formula the term containing the factor ordering, that plays a part only close 
to zero, is ignored. 

Because of the smallness of the sum J2k Cf under the barrier this term can be neglected 
in the expression for F e , and we suppose with a necessary accuracy 

F e (a) ~ -2a. (26) 

At that a the phase $ e has the form 

$ e = -/cm ( — 
\a 

and leaving only growing modes in solutions at a > a we obtain 




Figure 1: Form of the potential barrier without modes influence (dashed line) and with 
taking them in account (solid line). 



In that solution the amplitude of short waves catastrophically increases with growth of 
k but not under the exponential law. The latter is realized in the case of a slow variation 
of F e in eq. fl25|) . It is similar to the Rubakov-effect of catastrophic particle creation, see 
ref. g. 

A final value of an amplification of the amplitude of short waves is determined by a 
"time" of the WF stay under the barrier. Let us notice that we consider a self-consistent 
problem of amplification of short waves (an analogue of particle creation). It means that 
the energy of the increasing modes under the barrier controls the form of the barrier itself 
(see Fig. 1). 

Let us consider the equation for the action of the WF under the barrier (|13|) 
dS e V „ 2 Cl f a\ 2k 

Only a "potential" energy of increasing modes is present in the superpotential V, but 



here again it is not cancelled by their "kinetic" energy as it was realized in (|26|) . 

As it follows from V(ao, a±) = 0, an external boundary of the tunneling region is 
determined by the equation 

4 \a J 

The final value of amplification of the amplitude of the modes is defined by substitution 



of a± in eq. 

Further, the Euclidean action S e (ai) has a form 



Subject to k ^> 1 we have S e = —a\. Hence, a probability of tunneling of the WF through 
the barrier is given as 

w = exp(-2|S e |) = exp(-2a?) (30) 

at w -C 1. The formula 

exp(— 2a?) 

w = 

1 + exp(-2af) 

is more exact in the case w < 1. 

5 Conclusions 

The self-consistent problem of the evolution of the closed Friedmann Universe with an 
amplification of wave perturbations of a scalar field under a barrier is considered. Within 
the framework of a semiclassical approximation of the Wheeler-DeWitt equation it is 
shown that the increase of the amplitude of short waves under the barrier controls the 
barrier shape, i.e. leads to both a reduction of the barrier and fast output of WF from 
under the barrier. In this sense the problem is self-consistent. The process of tunneling 
and the increase of the amplitudes of short waves is defined here by the assignment of their 
initial spectrum in the inner classically allowed region. If it is required after tunneling in 
the Universe to receive a spectrum close to the flat one, the initial spectrum should be 
almost exponentially suppressed in the area k — > oo. It is qualitatively equivalent to the 
presence of an effective temperature in the inner classically allowed region. 

A more detailed discussion of these problems obviously requires a consideration of the 
quantum nature of the problem for short-wave perturbations of a scalar field. 
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